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(57) Abstract 



The invention relates to a method and a system for synthesizing a prescribed intensity pattern based on phase contrast imag'ng that 
is not based on the assumption of prior an methods that the pahase shift ^ is less than I radian. An improved method based on a simple 
imaging operation with a simple one-to-one mapping between resolution elements of a spatial phase modulator and resolution elements of 
the generated intensity pattern is provided. According to the invention a method is provided for synthesizing an intensity pattern with lou 
toss of electromagnetic energy, comprising spatial modulation of electromagnetic radiation with a spatial phase mask for modulation of the 
phase of the incident eletromagnetic radiation by phasor values of individual resolution elements of the spatial phase mask, each phasor 
value being detemiined in such a way that the values of the Fourier nansformed phason attain predetemiined values for predetermined 
spatial frequencies, and the phasor value of the specihc resolution element of the spatial phase mask corresponds to a distinct mirnMi> level 
of the image of the resolution element in the mteasity panem. and a spatial phase hiter for phase shifting of a pan of the electromagniic 
mdiatton. in combination uith an imagmg system for generation of the mtensii) panem h> interference in the image plane of the imat:ing 
\\%tem between the pan of the electromagnetic raidation that has been phase shifted b> the phase hIter and the remaining pan of ihe 
electromagnetic radiaiion. 
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Phase Contrast Imaging 
FIELD OP THE INVENTION 



PCT/DK96/00190 



The invention relates to a method and a system for synthesiz- 
ing a prescribed intensity pattern based on phase contrast 
imaging . 

BACKGROUND OF THE INVENTION 

It is well knovm to form an image on an illiominated surface 
of a body by absorption or blocking of energy of an illumi- 
nating beam. For example in an overhead projector, an over- 
head transparent absorbs or blocks part of the light beam of 
the projector whereby a large image of an overhead is formed 
on a screen. However, this results in a loss of light inten- 
sity as part of the emitted light from an image forming 
system is reflected or eOssorbed. 

To avoid loss of energy causing, e.g. loss of light intensity 
of the synthesized intensity pattern, power dissipation 
generating heat in components of the system, etc., methods 
and systems have been developed wherein the phase of a light 
beam is modulated instead of the amplitude or intensity of 
the light beam, as modulation of the phase of the light beam 
do not lead to loss of energy. The phase modulation is fol- 
lowed by a conversion of the phase modulation into an ampli- 
tude or intensity modulation. 

A diffractive optical element, such as a holographic optical 
element, may be used to generate a phase modulation. Then, 
the resulting intensity modulation at each point of a picture 
formed by conversion of the phase modulation into intensity 
modulation will depend upon the phase modulation values at 
each point of the diffractive optical element as the light 
intensity at each point of the picture is formed by a coher- 
ent superposition of light received from the entire surface 
of the diffractive optical element. Diffractive optical 
elements are rather complex to design for synthesis of a 
prescribed intensity pattern. 
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Imaging methods and systems may also be used in connection 
with phase modulation. These methods and systems are charac- 
terized by the fact that the intensity of a point of a pic- 
ture formed by conversion of phase modulation into intensity 
5 modulation will depend upon the phase modulation value of one 
point of the phase modulator only as this point is imaged 
onto the picture point in question by the imaging system* 
This one-to-one relationship makes the design of phase modu- 
lators in these systems simple. Methods and systems of this 
10 kind are named phase contrast imaging methods and systems. 

Phase contrast imaging methods were originally developed 
within the field of microscopy. Many objects of interest in 
microscopy are largely transparent, thus cd^sorbing little or 
no light. When light passes through such an object, the 
15 predominant effect is the generation of a spatially varying 
phase shift which can not be seen by a human as the eye of a 
human responds to light intensity and colour and does not 
respond to the phase of light. 

In 1935, Fritz Zemik proposed a phase contrast technique 
20 which rests on spatial -filtering principles and has the 

advantage that the observed intensity is linearly related to 
the phase shift introduced by the object. 

Suppose that a transparent object with amplitude 
transmittance 

25 t(x,y) = exp[j0fx,y;i (1) 

is coherently illuminated in an image -forming system. For 
simplicity, a magnification of unity is assumed and the 
finite extent of the exit and entrance pupils of the system 
is neglected. Further, a necessary condition to achieve 
30 linearity between phase shift and intensity is that the phase 
shift 4> be less than 1 radian, in which case the amplitude 
transmittance can be approximated by 
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t(x,y) =14 j4>(x,y) (2) 

The terms of order <t>^ and higher are neglected in this appro- 
ximation. It is seen that the first term of (2) leads to a 
strong wave con?)onent that passes through the sample without 
5 change, while the second term generates weaker diffracted 
light that is deflected away from the axis of the system. 

The image produced by a conventional microscope can be writ- 
ten 



I - |1 + j0|2 « 1 (3j 

10 where the term <t>^ has been approximated by zero. It is seen 
that the diffracted light is not observable because it is in 
phase quadrature with the strong background. As Zemik recog- 
nized that the background is brought to a focus on- axis in 
the focal plane while the diffracted light - containing 

15 higher spatial frequencies - is spread away from the focal 
point, he proposed that a phase - changing plate be inserted in 
the focal plane to modify the phase relation between focused 
and diffracted light. 



20 



The phase -changing plate can consist of a glass substrate on 
which a small transparent dielectric dot has been coated. The 
dot is placed at the center of the focal plane and has a 
thickness and index of refraction such that it retards the 
phase of the focused light by either 7r/2 radians or 37r/2 
radians relative to the phase retardation of the diffracted 
25 light. In the former case the intensity in the image plane 
becomes 

I = lexplj (7r/2)+j0|2 = |j(i+0)|2 -1 + 20 (4) 
while in the latter case 



I - jexplj (37r/2)+j0|2 « |j(-l*0)|2 •1-20 



(5) 
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Thus, the image intensity has become linearly related to the 
phase shift 0. When the phase of the background is retarded 
by 7r/2, the result is known as positive phase contrast, while 
a 3ir/2 retardation is said to yield negative phase contrast • 



5 It is seen that the method described above leads to a phase 
contrast imaging method that provides a small phase signal 
that is superimposed on a large DC- component . This leads to 
an important disadvantage of the method because, typically, 
it will be necessary to attenuate the DC- component to enhance 
10 the information contained in the phase modulated signal- 

However, the attenuation of the DC- component leads to loss of 
energy. This kind of filtering is usually denoted Dark Field 
Filtering. 

It is another disadvantage of the phase contrast imaging 
15 method described ahove that it is based on the assximption 
that the phase shift 0 is less than 1 radian which is very 
often not fulfilled in practical real -life applications. 
However, the theory is still applied to such applications, 
disregarding the fact that the basic assumption is not ful- 
20 filled, and this leads to non-optimized technical solutions. 

In EP 0 657 760 a phase contrast imaging system is disclosed 
in which an image simulation and projection system is based 
on the Texas Instrument flexure beam digital mirror device 
(DMD) . The flexure beam DMD is used for analog phase 

25 modulation of reflected light and the phase modulation is 

converted to amplitude modulation utilizing a phase contrast 
imaging method. The flexure beam DMD provides a flicker- free 
modulated wave and accordingly, optical image sensor 
synchronization is not needed. The system disclosed operates 

30 according to the Zernike method and, thus, includes the 
corresponding disadvantages described above. 



Another example of a phase contrast imaging system is 
disclosed in GB 2 199 716, wherein an optical guide -beam 
projector for a missile guidance system is disclosed that 
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provides a spatially intensity modulated guide-beam. A spa- 
tial phase modulator is used to generate the guide -beam. The 
phase encoding of the spatial phase modulator constitutes a 
periodic square-wave modulation (50% duty cycle) of two phase 
5 values 0 and 7r/2. The phase modulation is converted into an 
amplitude modulation by Fourier transforming lenses and a 
phase plate providing a phase shift of the background signal 
by 7r/2. A method for synthesizing the specific intensity 
pattern of the optical guide -beam based on phase contrast 
10 imaging is not disclosed in this document. 

A similar example of a phase contrast imaging system is 
disclosed in "Array illuminator based on phase contrast". 
Applied Optics Vol. 21, No. 14, pp. 2915-2921 (1988). A 
method is disclosed of converting a wide beam of uniform 

15 intensity into an array of bright spots without losses. The 
input spatial phase mask constitutes a periodic array of 
phase dots with the phase value tt, the remaining area of the 
phase mask having the phase value 0. The phase modulation is 
converted into an amplitude modulation by Fourier 

20 transforming lenses and a phase plate providing a phase shift 
of the background signal by tt. The method is limited to the 
implementation of periodic array configurations with the 
binary phase values 0 and tt. 

It is well-known to use so-called "radiation focusators", 
25 i.e. computer generated holographic optical elements, for 

spatial phase modulation of a light beam, e.g as disclosed in 
Special Issue on Computer Optics in the USSR, Optics and 
Lasers in Engineering, Vol. 15, no. 5 1991. However, such 
elements are complicated to synthesize. Typically, they are 
30 synthesized in such a way that the desired image is formed in 
the Fresnel region or the Frauenhofer region. Thus, the 
intensity of a resolution element in the generated image is a 
function of several, typically all, phase values of the 
resolution elements of the holographic optical element, 
35 Obviously, this complicates the design of a general purpose 
holographic optical element and advanced, very time consuming 
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algorithms have to be applied. Further, the complicated 
design of the holographic optical elements renders It almost 
Impossible to Implement dynamically changeable spatial phase 
modulators with such elements. 

5 It is a further disadvantage of holographic optical elements 
that a carrier frequency is needed to separate diffracted 
light from non-dif f racted light resulting in an off -axis 
system geometry and a need for a dif f ractive medium that can 
support these high frequency terms. 

10 SUMMARY OF THE INVENTION 

It is an object of the present invention to provide an appar- 
atus of the cibove kind which apparatus is robust, compact, 
simple to design and relatively cheap to manufacture. 

It is another object of the present invention to provide an 
15 improved method and apparatus for phase contrast imaging that 
take all terms of the Taylor's series: 

t(x,y) ^e^<f>^-'y) =i*i*(x,y) - » " -j » \ »(^^y) . . 

into account and, thus, is not based on the assumption that 
the phase shift 4> is less than l radian. It is important to 
20 note that each term of the Taylor's series contribute to the 
DC-value of the function t(x,y). This fact is not recognized 
in the art as the DC-value has until now been believed to be 
represented by numeral 1 in equation (2) . 

It is still another object of the present invention to pro- 
25 vide an improved method and apparatus for phase contrast 

imaging without the need of attenuating the DC- component of 
the signal to enhance the information contained in the phase 
modulated signal. 
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It is yet another object of the present invention to provide 
an in?>roved method based on a simple imaging operation with a 
simple one-to-one mapping between resolution elements of a 
spatial phase modulator and resolution elements of the gener- 
> ated intensity pattern. 



According to the invention a method is provided for syn- 
thesizing an intensity pattern with low loss of electromag- 
netic energy, comprising spatial modulation of electromag- 
netic radiation with a spatial phase mask for modulation of 
10 the phase of the incident electromagnetic radiation by phasor 
values of individual resolution elements of the spatial phase 
mask, each phasor value being determined in such a way that 

1) the values of the Fourier transformed phasors attains 
predetermined values for predetermined spatial fre- 
15 guencies, and 



) the phasor value of a specific resolution element of 
the spatial phase mask corresponds to a distinct 
intensity level of the image of the resolution 
element in the intensity pattern, 



and a spatial phase filter for phase shifting of a part of 
the electromagnetic radiation, in combination with an imaging 
system for generation of the intensity pattern by interfer- 
ence in the image plane of the imaging system between the 
part of the electromagnetic radiation that has been phase 
shifted by the phase filter and the remaining part of the 
electromagnetic radiation. 



Although, the present method is related to encoding of spa- 
tial phase masks in two spatial dimensions (planar encoding) , 
the principles of the method may be utilized for phase encod- 
ing in one to three spatial dimensions and/or in the temporal 
dimension. 
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The electromagnetic radiation may be of any frequency range 
of the electromagnetic spectirum, i.e. the gaimma frequency 
range, the ultraviolet range, the visible range, the infrared 
range, the far infrared range, the X-ray range, the microwave 
5 range, the HF (high frequency) range, etc. The present method 
is also appliccJDle to particle radiation, such as electron 
radiation, neutron radiation, etc. 



Preferably, the electromagnetic radiation is monochromatic or 
quasi -monochromatic so that the energy of the electromagnetic 

10 radiation is concentrated in a narrow frequency bandwidth. As 
the intensity pattern is synthesized by interference of two 
electromagnetic waves emitted from a common source of elec- 
tromagnetic radiation but the phases of which have been 
changed differently, it is required that the frequency range 

15 of the emitted electromagnetic radiation is sufficiently 

narrow to ensure that the two waves of electromagnetic radi- 
ation are coherent so that their superposition generates the 
desired intensity pattern. If the frequency range is too 
broad, the two waves will be incoherent and the phase infor- 

20 mation will be lost as superposition of non- coherent waves 
results in a summation of the intensities of the two waves. 
It is required that the difference between individual delays 
of electromagnetic radiation to be superpositioned is less 
than the wavelength of the radiation. This is a relaxed 

25 requirement that allows the electromagnetic radiation to be 
relatively broad-banded. For example in the visible range a 
Xe-lamp or a Hg-lamp can be used as a light source in a 
system according to the present invention with the advantage 
compared to a laser light source that the speckle noise is 

30 reduced. The requirements of the spatial coherence of the 
electromagnetic radiation depend upon the space bandwith 
product of the corresponding system and how close the 
required system performance is to the theoretically obtain- 
able performance of the system. 
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Preferably, the electromagnetic radiation is generated by a 
coherent source of electromagnetic radiation, such as a 
laser, a maser, a phase-locked laser diode array, etc. How- 
ever a high pressure arc lamp, such as a Hg lamp, a Xe lamp, 
5 etc, may also be used and even an incandescent lamp may be 
used as a source of electromagnetic radiation in a low per- 
formance system. 

A spatial phase mask is a component that changes the phase of 
an electromagnetic wave incident upon it. The spatial phase 
mask may transmit or reflect the incident electromagnetic 
wave. Typically, the spatial phase mask is divided into a 
number of resolution elements each of which modulates the 
incident electromagnetic wave by changing its phase by a 
specific predetermined value. The predetermined values are 
assigned to each resolution element in different ways depend- 
ing upon the technology applied in the component. For example 
in spatial light modulators, each resolution element may be 
addressed either optically or electrically. The electrical 
addressing technique resembles the addressing technique of 
solid-state memories in that each resolution element can be 
addressed through electronic circuitry to receive a control 
signal corresponding to the phase change to be generated by 
the addressed resolution element. The optical addressing 
technique addresses each resolution element by pointing a 
light beam on it, the intensity of the light beam correspon- 
ding to the phase change to be generated by the resolution 
element illuminated by the light beam. 

Spatial phase masks may be realized utilizing fixed phase 
masks, devices comprising liquid crystals and being based on 
30 liquid crystal display technology, dynamic mirror devices, 
digital micromirror arrays, deformable mirror devices, mem- 
brane spatial light modulators, laser diode arrays (inte- 
grated light source and phase modulator), smart pixel arrays, 
etc. 



10 



15 



20 



25 
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A spatial phase filter is typically a fixed phase mask, such 
as an optically flat glass plate coated with a dielectric 
layer at specific positions of the glass plate. However, the 
spatial phase masks mentioned in the previous section may 
5 also be used for spatial phase filters. 

» 

The imaging system maps the phase modulating resolution 
elements of the spatial phase mask on the target surface of 
the synthesized intensity pattern. It may comprise a 4f-lens 
configuration (two Fourier transforming lenses utilizing 
10 transmission of light or one Fourier transforming lens util- 
izing reflection of light) or a single imaging lens. However, 
any optical imaging system providing a filtering plane for 
the spatial phase filter may be applied in a phase contrast 
imaging system. 

15 In the method according to the present invention, the syn- 
thesized intensity pattern is generated by superposition of 
two electromagnetic waves in the image plane of the imaging 
system. The spatial phase mask changes the phase values of an 
electromagnetic wave incident upon it and the imaging system 

20 directs the electromagnetic wave with changed phases 

reflected from or transmitted through the spatial phase mask 
towards the spatial phase filter. The phase filter phase 
shifts a part of the electromagnetic radiation and the imag- 
ing system is adapted to superimpose in the image plane the 

25 phase shifted part of the electromagnetic radiation with the 
part of the electromagnetic radiation that is not phase 
shifted by the spatial phase filter. 

According to a preferred embodiment of the invention, the 
spatial phase mask is positioned at the front focal plane of 

30 a lens while the spatial phase filter is positioned in the 

back focal plane of the lens, whereby a first electromagnetic 
field at the phase mask is Fourier transformed by the lens 
into a second electromagnetic field at the phase filter. 
Thus, specific spatial frequencies of the first electromag- 

35 netic field will be transmitted through the spatial phase 
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filter at specific positions of the phase filter. For 
instance, the energy of the electromagnetic radiation at zero 
frequency (DC) is transmitted through the phase filter at the 
intersecting point of the Fourier plane and the optical axis 
5 of the lens also denoted the zero- order diffraction region* 

It is presently preferred that the spatial phase filter is 
adapted to phase shift the DC -part of the electromagnetic 
radiation and to leave the remaining part of the electromag- 
netic radiation unchanged or, alternatively, to leave the DC- 

10 part of the electromagnetic radiation unchanged and to phase 
shift the remaining part of the electromagnetic radiation. 
The last alternative is preferred when the energy level of 
the DC-part of the electromagnetic radiation is so high that 
the phase shifting part of the phase filter will be destroyed 

15 by it. For example in laser cutting, the DC level of the 
laser beam can be so high that a phase shifting dot posi- 
tioned at the intersecting point of the DC part of the laser 
beam at the phase filter would evaporate. It is also possible 
to block the electromagnetic radiation (no transmittance) in 

20 the zero-order diffraction region, however, the DC energy of 
the radiation is then lost. 

Below, an expression of the intensity of the synthesized 
intensity pattern as a function of the phasor values <t>(x,y) 
of the phase mask, when the DC-part of the electromagnetic 
25 radiation is phase shifted, is deduced. 

Electromagnetic radiation incident on the spatial phase mask 
can be described by a function A(x,y) , where A(x,y) is a 
complex number (amplitude and phase) of the incident field on 
the point (x,y) of the spatial phase mask. At the point 

30 {x,y) , the spatial phase mask modulates the phase of the 
incident radiation with a value 4>(x,y) so that the field 
after reflection by or transmission through the spatial phase 
mask may be described by the function A(x,y) ♦ e^*^*^'^^, 
^iit><x,y) being the phasor value of the point (x,y) of the 

35 spatial phase mask. As A(x,y) preferably is a constant value 
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over the entire surface of the spatial phase mask, the term 
is left out of the following equations for simplicity. 

The expression of the electromagnetic radiation incident on 
the spatial phase filter may now be separated into an AC- term 
5 and a DC- term. If the DC- term of the field is denoted a, the 
AC-term of the field is given by the term e^^^^'^^ -a. As the 
spatial phase filter changes the phase of the DC-part of the 
electromagnetic radiation by 6, the intensity of the syn- 
thesized intensity pattern at the image plane of the imaging 
10 system is given by: 

I(x\yM = \e^^<^''y'> + a{ei^-l)|2 (7) 

wherein {x',y') is the coordinates of the image of the point 
(X/y) of the spatial phase mask formed by the imaging system 
in the image plane. 

15 It should be noted that the second term of the equation is a 
complex niunber that adds to the phasors e^^(^*y> of the spa- 
tial phase mask and may be interpreted as a contrast control 
parameter for the synthesized intensity pattern I(x',y'). 

According to a preferred embodiment of the invention, the 
20 average value of the phasors is adjusted in order to control 
the range of intensity levels. 

Instead of phase shifting the DC-part of the electromagnetic 
radiation, it is also possible to synthesize a prescribed 
intensity pattern by phase shifting other parts of the elec- 

25 tromagnetic radiation by adapting the spatial phase filter to 
phase shift electromagnetic radiation incident upon one or 
more arbitrary regions of the phase filter and leaving the 
phase of the remaining part of the electromagnetic radiation 
unchanged and then superimposing the two parts of the elec- 

30 tromagnetic radiation. The corresponding mathematics and the 
corresponding design procedures for the spatial phase mask 
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and spatial phase filter will of course be more conplicated 
than for the method described in the previous section. 

A simple exanple of phase shifting a part of the electromag- 
netic radiation of a spatial frequency different from the 
5 zero frequency is provided by moving the DC-part of the 

electromagnetic radiation to another spatial frequency in the 
Fourier plane (identical to the plane of the spatial phase 
filter) utilizing an optical component with an appropriate 
carrier frequency (i.e. a grating or a prism) or, preferably, 
10 encoding the function of a grating or a prism into the spa- 
tial phase mask, and adapting the spatial phase filter to 
change the phase of the electromagnetic radiation at this 
spatial frequency and to leave the phase of the remaining 
part of the electromagnetic radiation unchanged. 

According to another preferred embodiment of the invention, 
the phase mask is not positioned in the back focal plane of 
the lens but in the Fresnel region of the lens instead. In 
this case, the electromagnetic field at the phase filter will 
be given by a Fresnel transformation of the electromagnetic 
field at the spatial phase mask. This further con5>licates the 
mathematics and the design procedures, for example the term o 
in equation (7) has to be substituted by the value of the 
Fresnel transformation at the point (s) of phase changes of 
the phase filter. However, the Fresnel transformation may be 
calculated from a Fourier transformation by multiplication of 
the phasor values of the spatial phase mask by a quadratic 
phase factor followed by a Fourier transformation. 

It is an important aspect of the present invention that each 
intensity level of the synthesized intensity pattern for each 
30 resolution element may be generated by at least two differen*-. 
phasor values of a resolution element of the spatial phase 
mask. 



15 



20 



25 
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For example, when the spatial phase filter phase shifts the DC- 
part of the electromagnetic radiation, it will be shown later 
that, advantageously, the average a of the phasors of the 
resolution elements of the phase mask should be eejual to M 
5 and the value of the phase shift 0 should be equal to tt. In 
this case, the intensity of the synthesized image pattern at 
the image (x',y') of the resolution element {x,y) will be 
given by: 

l{x',yn = 2(1 - cob4>(x' ,y')) (8) 

10 It is seen that con^lex conjugate phasors (values of 4> of 

opposite sign) result in identical intensity levels I(x',y')- 
It can be shown that for any value of the modulus of the 
average of the phasors |Qf|, two phasors exist that will 
generate identical intensity levels of the synthesized inten- 

15 sity pattern. 

Further, if the spatial phase filter phase shifts parts of 
the electromagnetic radiation different from the DC-part, the 
phasor value that generates a specific intensity level will 
depend on the position of the resolution element in question, 

20 i.e. the phasor value and the position of the resolution 

element with that phasor value together define the intensity 
level at the image of the resolution element in the syn- 
thesized intensity pattern. Still, it is true that for each 
resolution element of the spatial phase mask, each intensity 

25 level of the synthesized intensity pattern may be represented 
by one of two different phasors of complementary phase 
values . 

This freedom of being able to select, for each intensity 
level to be generated and for each resolution element of the 
30 spatial phase mask, one of two phasors is used to control the 
phase of the 'Fourier transform of the phasors at specific 
spatial frequencies by selection of phasors with appropriate 
phase values to ensure two intervals of biunique functional 
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dependence between phasor values and corresponding intensity 
values • 

This freedom of choice of phasors may be utilized to select 
phasors of neighbouring resolution elements of the spatial 
5 phase mask with a maocimum difference between them, thereby 
generating an electromagnetic radiation emitted from the 
phase mask with a maximum content of high spatial frequencies 
which will generate a good separation of the DC part of the 
electromagnetic radiation from its AC part. However, any 
10 other strategy of selecting between two possible phasor 

values of each resolution element may be chosen to generate a 
desired spatial frequency content of the electromagnetic 
radiation. 



Preferably, the phase of the Fourier transform of the phasors 
15 at specific spatial frecpaencies is adjusted in order to 
control whether the relation between each phasor and the 
corresponding intensity level is a monotonic increasing or a 
monotonic decreasing function. 

Below, a set of different methods are described that are 
20 provided according to the present invention for adjustment of 
the modulus of the Fourier transform of the phasors at speci- 
fic spatial frequencies to attain a prescribed value. If 
convenient, the methods may be combined. 



According to one of the methods, the individual phasors of 
25 the resolution elements of the phase mask are adjusted by a 
constant value until the desired value of the modulus of the 
Fourier transform of the phasors at specific spatial fre- 
quencies is attained while maintaining prescribed relative 
intensity levels between intensities of resolution elements 
30 of the intensity pattern, i.e. iteratively. 

According to another method, the individual phasors of the 
resolution elements of the phase mask are adjusted utilizing 
histogrcun techniques known from image processing. A histogram 
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is a bar chart showing the number of resolution elements of 
the synthesized intensity pattern with a specific intensity 
value as a function of the intensity value. Any histogram 
technique, such as histogram equalization, adapting the 
5 histogram to a predetermined distribution, etc., may be used 
iteratively until the modulus of the Fourier transform of the 
phasors at specific spatial frequencies attain the prescribed 
value . 

According to yet another method, the phasor pattern of the 
10 phase mask is spatially scaled in order to adjust the modulus 
of the Fourier transform of the phasors at specific spatial 
frequencies . 

According to still another method, the modulus of the Fourier 
transform of the phasors at specific spatial frequencies is 
15 adjusted utilizing half tone coding techniques, such as 
raster techniques, area ratio modulation, spot diameter 
modulation, etc. 

It is seen from the description above that the intensity 
levels may differ from one synthesized intensity pattern to 

20 the next as a consequence of the adjustments of the modulus 
of the Fourier transform of the phasors at specific spatial 
frequencies. Thus, it is preferred to control the power of 
the radiation source in dependence of the intensity range of 
the intensity pattern so that a sequence of different inten- 

25 sity patterns show uniform intensity levels. 

According to a preferred embodiment of the invention, the 
shape of the phase filter is adapted to match the spatial 
frequency content of the phasors of the spatial phase mask, 
e.g. to optimize the desired separation of the part of the 
30 electromagnetic radiation to be phase filtered from the 
remaining part of the electromagnetic radiation. 
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It is within the scope of the present invention that the 
imaging system further comprises zooming means for variable 
scaling of the synthesized intensity pattern. The zooming of 
the imaging system may be dynamically controllable, e.g. in 
5 response to the scaling of the pattern of phasor values of 
the phase mask. 

According to the present invention, the power of the radi- 
ation source may be controllable in response to the spatial 
scaling of the pattern in the phase mask and/or the zooming 
10 of the focusing system. 

in order to provide a compact and integrated system according 
to the present invention, the optical function of a Fourier- 
transforming lens is encoded into the phasors of the spatial 
phase mask. The Fourier transforming lens may be refractively 
15 or diffractively encoded into the phase mask. 

Similarly, the optical function of an output lens may be 
encoded into the phase filter either refractively or 
diffractively. 

Further, a condensation may be encoded into the phasor values 
20 of the spatial phase mask so that part of the electromagnetic 
radiation modulated by the phase mask has a substantially 
flat intensity profile in the image plane. Without this 
compensation, part of the electromagnetic radiation modulated 
by the phase mask will have a flat profile with perturbations 
25 resulting from the phase filtering superpositioned upon it. 
This may cause "ringings" (oscillations) at the edges of the 
synthesized intensity pattern. 

According to another preferred embodiment of the invention, 
the source of electromagnetic radiation comprises one or more 
30 light sources of different wavelengths corresponding to three 
different colours, such as red, green and blue, for gener- 
ation of intensity patterns of arbitrary colours. Further, 
several independent systems each one illuminated by its own 
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wavelength can be combined into a single multi* wavelength 
system. 

BRIEF DESCRIPTION OF THE DRAWINGS 



Fig. 1 shows a 4f optical system for phase contrast imaging, 

5 Fig. 2 shows a 2f optical system for phase contrast imaging, 

Fig. 3 shows a If optical system for phase contrast imaging, 

Fig. 4 shows (A) off -axis read- out of reflective SLM and (B) 
on -axis read -out of reflective SLM. 

Fig. 5 shows schematically an example of a prescribed inten- 

10 sity pattern in ID. 



Fig. 6 shows schematically the resulting phase encoding 
corresponding to Fig. 5. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 



Fig, 1 shows a 4f phase contrast imaging system (1) . A laser 
15 (2) emits a light beam which is expanded by a beam expander 
(3) into a plane light wave of uniform intensity and directs 
it towards a spatial phase mask (4) . The light beam is trans- 
mitted through the spatial phase mask (4) and a Fourier 
transforming lens (5) . The spatial phase mask is positioned 
20 in the front focal plane of the lens (5) and a spatial phase 
filter (6) is positioned in the back focal plane of the lens 
(5) that is also the front focal plane of a lens (7) . The 
Fourier transforming lenses (5, 7) need not have identical 
focal lengths. Different focal lengths lead to a magnifica- 
25 tion ratio different from one. The phase filter (6) phase 

shifts the zero order diffraction part (8) of the light phase 
modulated by the spatial phase mask (4). The synthesized 
intensity pattern is generated in the back focal plane (9) of 
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the lens (7) and a dynamic focusing system (10) images the 
synthesized intensity pattern onto a focusing plane (11) . 

The optical system is controlled by a computer (12) . The 
coinputer (12) comprises interface means for addressing each 
5 of the resolution elements of the phase filter (4) and trans- 
mitting a phasor value to the addressed resolution element. 
Further, the computer (12) comprises laser control means for 
controlling the power of the laser (2) and imaging control 
means for controlling the focusing and the image ratio of the 

10 dynamic focusing system (10) . The computer (12) also com- 
prises input means, such as a keyboard, a diskette drive, an 
optical disc drive, a network interface, a modem, etc, for 
receiving an image pattern to be synthesized by the system 
(1) . From the received image pattern, the computer is adapted 

15 to calculate phasor values to be transmitted to the resol- 
ution elements of the phase mask, e.g. based on a histograun 
technique as described herein. Optionally, the phase shift of 
the phase filter (6) is adjustable and controllable by 
optional phase control means of the computer (12) which may 

20 be further adapted to adjust the phase shift, e.g. utilizing 
equation (18) . 

Fig. 2 shows a 2f phase contrast imaging system (20) . A laser 

(21) emits a light beaun which is expanded by a beam expander 

(22) into a plane light wave of uniform intensity and directs 
25 it towards a spatial phase mask (23) and a polarization beam 

splitter (24) and a quarter- wave plate (25) . The polarization 
beam splitter (24) and the quarter-wave plate (25) allows 
beam- splitting of light of a specific linear polarization 
without the power loss associated with conventional beami- 

30 splitters due to splitting of the beam in both directions of 
transmission through the beam-splitter. After transmission 
through the polarization beam splitter (24) and the quarter- 
wave plate (25) , the light beam is transmitted through a 
Fourier transforming lens (26) and is reflected from a spa- 

35 tial phase filter (27). The spatial phase mask (23) is posi- 
tioned in the front focal plane of the lens (26) and the 
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spatial phase filter (27) is positioned in the back focal 
plane of the lens (26) , The phase filter (27) phase shifts 
the zero order diffraction part (28) of the light that is 
phase encoded by the spatial phase mask (23) • The synthesized 
5 intensity pattern is generated in the back focal plane (29) 
of the lens (26) and a dynamic focusing system (30) images 
the synthesized intensity pattern onto a focusing plane (31) . 
As described for the system shown in Fig. 1, the system (20) 
is controlled by a computer (32) . 

Fig- 3 shows a If phase contrast imaging system (40) . A laser 

(41) emits a light beam which is expanded by a beam expander 

(42) into a plane light wave of uniform intensity and directs 
it towards a spatial phase mask (43) , The light beam is 
transmitted through the spatial phase mask (43) and an image 
forming lens (44) . A phase filter (45) positioned in the back 
focal plane of the lens (44) phase shifts the zero order 
diffraction part of the light phase encoded by the spatial 
phase mask (43) . The synthesized intensity pattern is gener- 
ated in the image plane (46) of the lens (44) and a dynamic 
focusing system (47) images the synthesized intensity pattern 
onto a focusing plane (48) . As described for the system shown 
in Fig. 1, the system (40) is controlled by a computer (49). 

Pig. 4 shows details of (A) an off -axis read-out of a reflec- 
tive phase mask (50) (or a spatial light modulator) and of 
25 (B) an on-axis read-out of a reflective phase mask (51) with 
a beam splitter (52) . Both configurations (A, B) may be 
utilized in the systems shown in Figs. 1-3. 



10 



15 



20 
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PHASE ENCODING FOR DC PHASE FILTERING 

In the following an example of encoding a spatial phase 
mask and a spatial phase filter will be given based on a 
5 system filtering in the DC-frequency range. The system 
chosen in this example is based on a 4-f lens 
configuration as shown in Fig. 1 and illuminated by 
electromagnetic radiation in the visible frequency 
domain, hereafter simply denoted as light radiation. 



Assuming that the illuminating light is monochromatic and 

has a substantially flat amplitude profile we obtain the 

following spatial amplitude distribution emitted from the 
spatial phase mask: 



where a{x, y) = exp(i^x, y)) represent the spatially encoded 
phasor values and AxAy is the area of the input phase 
20 modulating spatial light modulator. 

It turns out to be convenient to separate a{x, y) into two 

terms describing a spatially invariant DC-value, a, and a 

spatially varying AC-contribution Aa(x, y) . The DC-value 
23 can be found as: 



10 



15 





(10) 



Subsequently the AC-term is expressed by: 



30 



Aa(x. y) = exp{i(!H>:, y)) — ff exp(i(«x, y))dxdy (11) 

AxA\' ^•'j; 



wo 96/34307 



22 



PCT/DK96/00190 



The separation of a{x, y) into a spatially invariant DC- 
term and a spatially varying AC-term is an important 
point and will be used throughout the remaining part of 
5 this example, especially in the description of the 
spatial filtering procedure. 

The spatial filter utilized in this example is chosen as 
a circular phase contrast filter (different transverse 
10 shapes can also be used) centered around origo in the the 
spatial frequency domain, denoted by coordinates (f^, fy) : 

r(/,) = 1 + (exp(i^) - l)circ^-^| (12) 

15 where jf, = ^f^ + f, denotes radial spatial frequency and 
A/, describes the size of the circular (circ) phase 
filter. 



20 



In the spatial frequency domain (the filtering plane) the 
Fourier transformation (3) of the spatially modulated 
light radiation from the spatial phase mask is present. 
The filtering operation on the Fourier transformed light 
radiation performed by the spatial phase contrast filter 
can be expressed as a simple point-by-point 
25 multiplication procedure. Subsequently the spatially 
filtered light is inverse Fourier transformed ( TM by the 
second Fourier lens (Fourier transformation and reflected 
output coordinates) and the resulting spatial amplitude 
distribution in the image plane (with coordinates (x', y') ) 
30 can accordingly be written as: 
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o(x',y') = a(xVy') + (exp(i^) - l)3-|3(a(x, y))circ|^-^jj 

= a{x', y') + a(exp(i^) - l)rect^^ , j (13) 
= [exp(i^Kx', y')) + a(exp(i^) - l)]rect^-£ , -^j 

Within the illumination-region, (x', y') 6 , outlined by 

recti — ,—], one obtains: 
^Ax Ay J 



|o(x', y')f = 1 + 4p| sin[|J|q sin[|) - sin(^, - <^x', y') |]] (14 



Requiring that |o(x;, y;)f s 0 corresponding to complete 
darkness as the lowest intensity level in 
10 regions (x;, y^) e SR; implies : 



1 + 



^ sin(UH sin(|] - sin(^, - ^ |]] = 0 (15) 



where the abbreviation ^„ = ^x;, y^) has been used. 

IS 

The solutions to Eq. (15) are given by: 



20 



SX 

H = r . (e\ 



The requirement o < ja| < 1 implies that 
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e - n - - + pZ/ff P = 0,±1,±2, . . , 



leading to 



2 sin] 

5 



iij 



(18) 



where the +sign is for ^-values in the Interval: 



^ €]|;^ + P_2;r (19) 



10 and the -sign is for ^-values: 

The corresponding interval for - is: 

15 

^ ]f f[ <2i' 

Inserting the expression for |a| , one obtains the simple 
intensity expression: 

20 

' '^^ (22) 



where 



|o(x', y'i = 2|^1 T sin^<^5 - <*(x', y') + |j 



JJ exp(i ^x, y))dxdy = AxAy|a| exp(i ^- ) (23) 
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The phase-only transformations imply that energy is 
conserved: 

JJ |o(xV y')f cfx'dy' = jj \a{x, y)fdxciy = AxAy (24) 



A special case: 

The most convenient choice for a is: ^ " ^ (implying 

that ^ = ;r + Peven2^)' ^o that the output intensity can be 
10 described as: 



|o(x'.yf = 2[l - cos(((Kx', y'))] (25) 

In this case the phase-^intensity mapping is described by 
15 the intervals [O; ;r] -> [0;4] . 

By setting ^ = ^ one obtains the following requirement to 

the phase function ^x, y) : 

fj co^^x, y))dxdy = ^ 

^ (26) 

Jl sin(^x, y))dxdy = 0 

20 

Inserting the expression for \c(x', y')f in Eq. (24) yields: 

2 jj [l - cos(^x', y'))]:ix'dy' = AxAy (27) 



25 m accordance with the first of the integral expressions 
in Eq. (26) . 
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Encoding procedure: 

• A given intensity distribution (image) |o(x'/ y'^ is 
desired at the output side of the optical setup. 

5 • Pixellation of the image, that is generally 
represented in the greyscale range: [0;gmax], provides 
the relation: 
JJ |o{x', y')f dx'dy' = AxAy => ^ |o(i, j)f = Pi^(^y) • 



10 



The histogram for the desired image |o(i, j)f is 

adjusted (adj) within the greyscale range [0:gmax], so 
that the previous point is fulfilled: 

|o(i, j)f |o(i, . 



15 • The phase values can now be calculated as 

/ ...2 A 



^i/ j) = arccos 



1 - 



ginax 



20 



As before pixellation provides the relation; 

X sin((^i, j)) = 0 . 



•«3 



The previous point can now be fulfilled by complex 
conjugating half the input pixels having the same 
phase value in the phase histogram. 



25 • The phase conjugate phase flipping provides a valuable 
tool (an extra degree of freedom) for manipulating the 
spatial frequency content in order to optimize the 
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separation of low and high frequency terms at the 
filter plane. 

• The scheme is robust to constant phase errors across 
5 the input spatial phase modulator/ since Eq. (22) is a 

function of the difference: - j) f only. 

Furthermore, small variations in the individual pixel 
phase values do not introduce any detrimental effects 
because the average value a , is a result of a very 
10 large phasor sum. 

• If the desired intensity distribution is too small to 
include all energy, that is, the histogram is scaled 
to maximum and the left hand side of Eq. (24) is still 

15 smaller than the right hand side, then the input phase 

object can be scaled until Eq, (24) is fulfilled. In 
order to obtain a scale invariant output intensity 
level a dynamic focusing system is needed. Similarly, 
intensity invariance can be obtained by controlling 

20 the radiated power from the light source. 

Alternatively, one can ignore the residual background 
illumination and obtain intensity levels with a gain 
factor of 9" (bacJcground constant equal to 1") for 
narrow generally shaped line structures (e.g. Eq. 

25 (14)). 

Example 1: 

A very simple example illustrating the individual steps 
in the above procedure will be given below. To simplify 
30 the example it will be considered in one dimension only. 
The starting point for encoding the spatial phase mask in 
this example is based on the following parameters: 
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10 



15 



20 



a = 0.5 

^ = (28) 
gmax - 4 

Consider the pixellated 3-step function shown in Fig. 5 
to be synthesized in the image plane as an intensity 
5 distribution. From the above choices of parameters one 
obtains the simple relation between phase values in the 
spatial phase mask and the image intensity values: 

|o(i)f = 2[l - cos{M)] (29) 



To proceed from here it necessary to calculate the 
accumulated intensity S W-^^ image to be 

i 

synthesized. The accumulated intensity is easily 
calculated from an image histogram where the x-axis 
represents greylevel value and the y-axis represents the 
amount of pixels in the image at a given greylevel value. 

By use of a histogram ^ |o(i)f is simply found as the 

weighted sum of all greylevel values (x-axis) multiplied 
by their pixel counting (y-axis) . This describes, so to 
speak, the "weight" of the image. In this simple example 
histogram calculations are not needed since we only have 
3 greylevels with well-defined separations. 

The value for the accumulated intensity has to obey the 
25 equality: 



From Fig. 5 we obtain: 
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2 |o(i)r = 4pixels • 0 + 4pixels • (0.5max) + 6pixels • max = dmax 

i 

(31) 

5 So that the value for max can be estimated to be: 

7 

max = - (32) 
4 

The corresponding adjusted intensity levels, |o(i)|^^, are 

10 therefore: 7/4, 7/8 and 0. These values can now be 
utilized to calculate the phase values of the spatial 
phase mask from the relation: 



( 2|o(i)| 
^i) = arccos 1 - 



gmax 



= arccos 



(33) 



15 



where from we obtain the three phase values: 1.4 5 rad. 
0.97 rad. and 0 rad. 

The last step needed in order to encode the spatial phase 
mask is that the following equality is fulfilled: 



20 



2 sin{M) = 0 



(34) 



Since we have the choice to use complex conjugate phasor 
values (two phasors giving the same intensity level) many 
25 approaches can be taJcen from here. A simple approach is 
to flip every second phasor with its complex conjugate 
value as shown in Fig. 6. The final phase values used in 
the phase mask are accordingly: ± 1.46 rad. ± 0.97 rad. 
and 0 rad. 
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As the last step we can check whether the criteria: 
a = 1/2, is actually fulfilled with the chosen phasor 
encoding: 

o = — (4 exp(iO) + 2 exp(i0.97) + 3 exp(iL45) + 2 exp(-i0^7) + 3 exp(-iL45)) « 1 / 2 
14 

(35) 
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GENERAL PHASE CORRECTION PROCEDURE INTEGRATED WITH THE 
PHASE ENCODING 

In Eq. (14) we obtained an analytic relation between the 
phase values in the spatial phase mask and the resulting 
intensity distribution, within the region (x', y') e W : 



The analysis leading to the above relation was based on 

the assumption that |a| is a constant value within the 

domain. In other words, the following approximation was 
applied: 



However, for certain spatial filter parameters the 
lefthand side of this expression will not be a space 
invariant constant value throughout the whole 9?' -domain 
but will instead manifest slowly variations/oscillations. 
This will introduce small errors in the final 
superposition between the phase filtered DC-value and the 
direct propagated AC-signal, In order to circumvent this 
problem a technique is needed that can counteract the 
distortions by use of phase-only encoding in the 
components already present in the system. In what follows 
a procedure for integrating predistortion that 
counteracts the above mentioned distortions will be 
described that is purely based on modifying the phasor 
values in the spatial phase mask at the input side of the 
system. The method can also counteract other types of 
distortions inherent in a practical implementation of the 




(36) 



3(a(x, y))circ 




(37) 
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system. Furthermore, the method can be applied in systems 
filtering at other spatial frequencies than DC. 

Procedure: 

When encoding the input phase function it is helpful to 
have a reverse equation, expressing the input phase 
distribution as a function of an adjusted (electronic) 
image grey-level distribution, addressing the input 

spatial light modulator: 

-H^ = 1 + 4p(x'. y')| sin^llpCx', y')| sin[|) - sin^^^^^Cx', y') - y') + |]j 

(38) 

where it has been talcen into account that a(x', y') is not 
considered as a constant but manifests a smooth 
oscillating behaviour within the optical image domain. 
The maximum value of is denoted gmax . 

Now, one can derive a formula for the grey-level 
correction, AI,,„(x', y') . that one needs to apply in order 
to encode a phase function that compensates for the 
spatial variation of the average phase value a(x', y') : 

'^-"^''^'/^ £ 1 + ^(x; y')| sin[|][(5(x% y')! sin(|] - sin(«»,(x', y') - ^x'. y') * |]] 

J 2(j...(x', y) * AI.,.(x'. y)) \ 
^x-.y) = acosj^l — J 

(39) 

where the second relation has been derived from the first 
by setting a = 1 / 2 and $ = n . 

By inserting the second relation m the first expression 
one gets: 



BAD ORIGJNA^ 
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f \ 

A J.,,(x', y') = 



^— L--^ - 1 X„„(x', y') - ^=£?ll£ii^ tex', y')| - 
^2(a(x',y')| j 2|a(x', y')|V ' 2^ 



(40) 



5 This formula is however not directly useful because it is 
related to the histogram adjusted grey-level distribution 
denoted by J,j^ . 

One needs a formula that relates the above correction 
10 terra to the original input grey-level distribution J(x, y) 
that has not been modified by histogram adjustments. This 
is important since the effect of the grey-level 
corrections also have to be incorporated in the procedure 
of histogram adjustments. 

IS 

The histogram scaling gives: 



I(x, y) = I,,„(x, y) (41) 



20 where J^^ and I,i«,aax the maximum grey-level values 

occurring in the original and the adjusted electronic 
grey-level distributions respectively. 

Similarly, one can apply this relation to the intensity 
25 correction term AJ^^^ and obtain: 

I(x, y) = I(x, y) + Al(x, y) = (l,, Jx, y) + AI„„(x. y)) 

^sJa^max 

(42) 

resulting in: 

.^0 
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I(x, y) = -j=:^ fl(x, y) - gmax f|a(x, y)| - \1] (43) 

In order to have enough dynamic range in grey-levels for 
the correction term one can derive an inequality from the 
above relation by using the fact that I^^ < gmax : 



or 



2|amln| 


• gmax 


^ gjnax 







(45) 



Since the first term is the dominating term in the 

expression for the intensity correction it will in 

practice be sufficient just to have the much simpler 
corrections: 



^ 2la(x, y)| (4 6) 

S 2|a„i„| • gmax 
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PROPOSED APPLICATIONS: 

• Laser machining, marking, branding, trimming, hardening, 
5 scribing, labeling, welding and cutting on two- and 

three-dimensional surfaces especially by use of CO2 and 

Nd:YAG laser based systems. The main advantage is that 

energy is not absorbed in the system (thereby preventing 

damage of the optical hardware) and this nonabsorbed 

10 energy is instead utilized to increase the intensity 

level of the desired light distribution in the image 

plane. High power can be delivered to selected regions on 

a work piece simultanously . 

15 • Efficient and dynamic spot-array generators based on 
phase contrast imaging. In order to provide bias or 
holding beams for arrays of optoelectronic elements, such 
as bistable elements, photonic switches and smart pixels. 

20 • Generation of structured light (lossless) for machine 
vision applications. E.g. periodic and skew periodic mesh 
grid illumination that can be updated in parallel. 

• Photolithographic applications (laser 3D direct writing 
25 in parallel without the need for sequential scanning) . 

E.g. high power laser direct writing of waveguides in Ge- 
doped silica. 

• Spatial light intensity modulation in general by use of 
30 pure phase modulation (radiation focusators) . 



• Laser beam shaping (dynamic) . 
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• Highly efficient parallel image projection without the 
need for a laser scanning device. 

5 • Dynamic Infrared Scene Projection (DIRSP) . 

• Exposure device for grating and mask production. 

• LIDAR applications. 

10 

• Laserprinting in parallel. 

• Lasershow applications. 
15 • Atmosphere research. 
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1. A method for synthesizing an intensity pattern with low 
loss of electromagnetic energy, comprising spatial modulation 
of electromagnetic radiation with 
5 a spatial phase mask for modulation of the phase of the 

incident electromagnetic radiation by phasor values of indi- 
vidual resolution elements of the spatial phase mask, each 
phasor value being determined in such a way that 

1) the values of the Fourier transformed phasors attains 
10 predetermined values for predetermined spatial fre- 

quencies, and 



2) the phasor value of a specific resolution element of 
the spatial phase mask corresponds to a distinct 
intensity level of the image of the resolution 
15 element in the intensity pattern, and 



a spatial phase filter for phase shifting of a part of 
the electromagnetic radiation, in combination with 

an imaging system for generation of the intensity pattern 
by interference in the image plane of the imaging system 
20 between the part of the electromagnetic radiation that has 

been phase shifted by the phase filter and the remaining part 
of the electromagnetic radiation. 



2. A method according to claim 1, wherein the fact that for 
each resolution element of the spatial phase mask, each 

25 intensity level of the synthesized intensity pattern may be 
represented by one of two different phasors of complementary 
phase values is used to control the phase of the Fourier 
transform of the phasors at specific spatial frequencies by 
selection of phasors with appropriate phase values in such a 

30 way that the function defining the intensity values as a 

function of the phasor values comprises two regions of biuni- 
que functional dependence between phasor values and corre- 
sponding intensity values. 
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3. A method according to claim 1 or 2, wherein the modulus of 
the Fourier transform of the phasors at specific spatial 
frequencies is adjusted in order to control the range of 
intensity levels of the synthesized intensity pattern. 

5 4. A method according to any of claims 1-3, wherein the phase 
of the Fourier transform of the phasors at specific spatial 
frequencies is adjusted in order to control whether the 
relation between each phasor and the corresponding intensity 
level is a monotonic increasing or a monotonic decreasing 
10 function. 

5. A method according to any of claims 1-4, wherein the 
adjustment of the modulus of the Fourier transform of the 
phasors at specific spatial frequencies comprises at least 
one of the following measures: 

15 a) adjusting the individual phasors of the resolution 

elements of the phase mask maintaining prescribed relative 
intensity levels between intensities of resolution elements 
of the intensity pattern, 

b) adjusting the individual phasors of the resolution 
20 elements of the phase mask by histogram techniques, 

c) spatially scaling the phasor pattern of the phase 
mask, 

d) utilizing half tone coding techniques. 

6. A method according to any of claims 1-5, wherein the power 
25 of the radiation source is controlled in response to the 

intensity range of the intensity pattern. 
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7. A method according to any of claims 1-6, wherein each 
phasor of the phase mask is selected from the set of two 
phasors with complementary phase values in such a way that a 
specific spatial frequency distribution of the intensity of 

5 the electromagnetic radiation at the spatial filter is 
attained* 

8. A method according to claim 7, wherein the phase of phaso- 
rs of adjacent resolution elements and of substantially 
identical phasor values alternates between the two possible 

10 complementary phase values of phasors. 

9. A method according to claim 8, wherein the shape of the 
phase filter is adapted to match the spatial frequency con- 
tent of the phasors of the spatial phase mask. 

10. A method according to any of claims 1-9, wherein phasors 
15 of the phase mask and the phase shift of the phase filter are 

selected in such a way that the intensity level of at least 
one resolution element of the synthesized intensity pattern 
is zero. 

11. A method according to any of claims 1-10, wherein the 
20 modulus of the Fourier transform of the phasors at specific 

spatial frequencies is adjusted, by selection of appropriate 
phasors of the phase mask, substantially to its minimum 
values to obtain a maximum dynamic range of the phasors of 
the phase mask. 

25 12. A method according to any of claims 1-11, wherein the 
phase filter is centred around the zero-order diffraction 
region. 

13. A method according to claim 12, wherein at least one 
intensity level of the intensity pattern may be represented 
30 by one of two complex conjugate phasors. 
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14. A method according to claim 12, wherein at least one 
intensity level of the intensity pattern is represented by 
one of two phasors with phases the sum of which equals tt. 

15. A method according to any of claims 10-14, wherein the 
5 modulus of average value of the phasors and the phase shift 

introduced by the phase filter of the zero- order diffracted 
electromagnetic radiation are related to each other in accor- 
dance with the equations 



>°"=7r— T7 and 
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10 16. A method according to claim 15, wherein modulus of the 
average value of the phasors ranges from 0.1 to 0.9. 

17. A method according to claim 15 or 16, wherein modulus of 
the average value of the phasors ranges from 0.25 to 0.75. 

18. A method according to any of claims 15-17, wherein 

15 modulus of the average value of the phasors ranges from 0.4 
to 0.6. 

19. A method according to any of claims 15-18, wherein 
modulus of the average value of the phasors is approximately 
0.5. 

20 20. A method according to any of claims 1-19, wherein the 
phase shift of the phase filter ranges from 7r/4 to 77r/4. 

21. A method according to any of claims 1-20, wherein the 
phase shift of the phase filter ranges from ir/2 to 3Tr/2. 



22. A method according to any of claims 1-21, wherein the 
25 phase shift of the phase filter ranges from 37r/4 to 57r/4. 
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23 • A method according to any of claims 1-22, wherein the 
phase shift of the phase filter is approximately tt. 

24. A method according to any of claims 1-23, wherein the 
phase shift of the phase filter is variable, 

5 25. A method according to any of claims 1-24, wherein the 
imaging system comprises zooming means for variable scaling 
of the intensity pattern. 

26. A method according to claim 25, wherein the zooming of 
the imaging system is dynamically controllable. 

10 27. A method according to claim 25 or 26, wherein the zooming 
of the imaging system is controllable in dependence of the 
scaling of the phase mask. 

28. A method according to any of claims 25-28, wherein the 
power of the radiation source is controllable in response to 

15 the spatial scaling of the pattern in the phase mask and/or 
the zooming of the focusing system. 

29. A method according to any of claims 1-28, wherein the 
phase mask is a spatial light modulator. 

30. A method according to any of claims 1-29, wherein the 
20 optical function of a Fourier- transforming lens is encoded 

into the phasors of the phase mask. 

31. A method according to any of claims 1-30, wherein the 
optical function of an output lens is encoded into the phase 
filter. 



25 32. A method according to any of claims 1-31, wherein part of 
the electromagnetic radiation modulated by the phase mask has 
a substantially flat intensity profile in the image plane. 
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33. A method according to any of claims 1-32, wherein the 
source of electromagnetic radiation comprises one or more 
light sources of different wavelengths corresponding to three 
different colours, such as red, green and blue, for gener- 

5 ation of intensity pattern of arbitrary colours. 

34. A phase contrast imaging system for synthesizing an 
intensity pattern with low loss of electromagnetic energy, 
comprising 

at least one source of electromagnetic radiation for emission 
10 of electromagnetic radiation, 

a spatial phase mask for modulation of the phase of the 
incident electromagnetic radiation by phasor values of indi- 
vidual resolution elements of the spatial phase mask, wherein 
each phasor value has been determined in such a way that 

15 1) the values of the Fourier transformed phasors attains 

predetermined values for predetermined spatial fre- 
quencies, and 



2) the phasor value of a specific resolution element of 
the spatial phase mask corresponds to a distinct 
20 intensity level of the image of the resolution 

element in the intensity pattern, and 



a spatial phase filter for phase shifting of a part of the 
electromagnetic radiation, and 

an imaging system for generation of the intensity pattern by 
25 interference in the image plane of the imaging system between 
the part of the electromagnetic radiation that has been phase 
shifted by the phase filter and the remaining part of the 
electromagnetic radiation. 
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35. A system according to claim 34, wherein the imaging 
system con^rises a first and a second Fourier transforming 
lens I the spatial phase mask being positioned in the front 
focal plane of the first lens, the spatial phase filter being 

5 positioned at the back focal plane of the first lens, and the 
second lens being positioned so that its front focal plane is 
positioned at the position of the back focal plane of the 
first lens. 

36. A system according to claim 34 or 35, wherein the imaging 
10 system comprises one Fourier transforming lens the spatial 

phase filter being positioned at the back focal plane of the 
lens . 

37. A system according to any of claims 34-36, wherein the 
imaging system comprises one imaging lens the spatial phase 

15 filter being positioned in the back focal plane of the lens. 

38. A system according to claim 34-37, further comprising a 
polarizing beam splitter and a quarter wave plate and/or a 
phase filter reflecting electromagnetic radiation incident 
upon it. 

20 39. A system according to any of claims 34-38, wherein the 
spatial phase filter changes the phase of the DC-part of the 
electromagnetic waves and leaves the phase of the remaining 
part of the radiation unchanged. 

40. A system according to any of claims 34-38, wherein the 
25 spatial phase filter do not change the phase of the DC-part 

of the electromagnetic waves and changes the phase of the 
remaining part of the radiation. 

41. A system according to any of claims 34-38, wherein the 
spatial phase filter blocks the DC-part of the electromag- 

30 netic waves and leaves the remaining part of the radiation 
unchanged. 
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42. A system according to any of claims 34-41, wherein the 
source of electromagnetic radiation is a Laser. 
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